The Nuclear Level Density and the Determination of 
Thermonuclear Rates for Astrophysics 



T. RAUSCHER\ F.-K. THIELEMANN\ K.-L. KRATZ^ 

^ Institut fiir theoretische Physik, Universitdt Basel, Basel, Switzerland 

^Institutf. Kemchemie, Universitdt Mainz, Mainz, Germany 



ABSTRACT. The prediction of cross sections for nuclei far off stability is crucial in the field of 
nuclear astrophysics. In recent calculations the nuclear level density — as an important ingredient to 
the statistical model (Hauser-Feshbach) - has shown the highest uncertainties. We present a global 
parametrization of nuclear level densities based on the back-shifted Fermi-Gas formalism. Employment 
of an energy-dependent level density parameter a and microscopic corrections from a recent FRDM 
mass formula by MoUer et al. leads to a highly improved fit of level densities at the neutron-separation 
energy in the mass range 20 < A < 245. The importance of using proper microscopic corrections from 
mass formulae is emphasized. The resulting level description is well suited for astrophysical applications. 

The level density can also provide clues to the applicability of the statistical model which is only 
correct for a high density of excited states. Using the above description one can derive a "map" for the 
applicability of the model for reactions of stable and unstable nuclei with neutral and charged particles. 

1. Introduction 

Explosive nuclear burning in astrophysical environments produces unstable nuclei, which 
again can be targets for subsequent reactions. In addition, it involves a very large number 
of stable nuclei, which are not fully explored by experiments. Thus, it is necessary 
to be able to predict reaction cross sections and thermonuclear rates with the aid of 
theoretical models. Explosive burning in supernovae involves in general intermediate 
mass and heavy nuclei. Due to a large nucleon number they have intrinsically a high 
density of excited states. A high level density in the compound nucleus at the appropriate 
excitation energy allows to make use of the statistical model approach for compound 
nuclear reactions (e.g. Hauser & Feshbach 1952, Mahaux & Weidenmiiller 1979), which 
averages over resonances. 

As the capture of an alpha particle leads usually to larger Q-values than neutron 
or proton captures, the compound nucleus is created at a higher excitation energy and 
especially in the case of alpha-captures it is often even possible to apply the Hauser- 
Feshbach formalism for light nuclei. Another advantage of alpha-captures is that the 
capture Q-values vary little with the N/Z-ratio of a nucleus, for nuclei with Z<50. For 
Z>50, entering the regime of natural alpha-decay, very small alpha-capture Q-values 
can be encountered for proton-rich nuclei. Such nuclei on the other hand do not play a 
significant role in astrophysical environments, maybe with exception of the p-process. 
This means that in the case of alpha-captures the requirement of large level densities 



at the bombarding energy is equally well fulfilled at stability as for unstable nuc;lei. 
Opposite to the behavior for alpha-induced reactions, the reaction Q-values for proton 
or neutron captures vary strongly with the N/Z-ratio, leading eventually to vanishing 
Q-values at the proton or neutron drip line. For small Q-values the compound nucleus 
is created at low excitation energies and also for intermediate nuclei the level density 
can be quite small. Therefore, it is not advisable to apply the statistical model approach 
close to the proton or neutron drip lines for intermediate nuclei. For neutron captures 
close to the neutron drip line in r-process applications it might be still permissable for 
heavy and often deformed nuclei, which have a high level density already at very low 
excitation energies. 

In astrophysical applications usually different points are emphasized than for in- 
vestigations in pure nuclear physics. Many of the latter in this well established field 
were focused on specific reactions, where all or most "ingredients", like optical poten- 
tials for particle and alpha transmission coefficients, level densities, resonance energies 
and widths of giant resonances to be implementated in predicting El and Ml gamma- 
transitions, were deduced from experiments. This of course, as long as the statistical 
model prerequisits are met, will produce highly accurate cross sections. For the majority 
of nuclei in astrophysical applications such information is not available. The real chal- 
lenge is thus not the well established statistical model, but rather to provide all these 
necessary ingredients in as reliable a way as possible, also for nuclei where none of such 
informations are available. In addition, these approaches should be on a similar level as 
e.g. mass models, where the investigation of hundreds or thousands of nuclei is possible 
with manageable computational effort, which is not always the case for fully microscopic 
calculations. 

The major uncertainty in all existing calculations stems from the prediction of nu- 
clear level densities (Truran et al. 1966; Holmes ct al. 1976: Cowan, Thielemann & 
Truran 1991, and extended references therein), which in earlier calculations gave uncer- 
tainties even beyond a factor of 10 at the neutron separation energy (Gilbert & Cameron 
1965), about a factor of 8 (Woosley et al. 1978), and a factor of 5 even in the most re- 
cent calculations (e.g. Thielemann, Arnould & Truran 1987; see also Fig. 3.16 in Cowan, 
Thielemann & Truran 1991). 

2. Thermonuclear Rates from Statistical Model Calculations 

A high level density in the compound nucleus permits to use averaged transmission 
coefficients T, which do not reflect a resonance behavior, but rather describe absorption 
via an imaginary part in the (optical) nucleon-nucleus potential (for details see Mahaux 
& Weidenmiiller 1979). This leads to the well known expression 
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for the reaction i^{j,o)m'^ from the target state to the exited state m'^ of the final 
nucleus, with center of mass energy Ey and reduced mass /iy . J denotes the spin, E the 
excitation energy, and tt the parity of excited states. When these properties are used 
without subscripts they describe the compound nucleus, subscripts refer to states of 
the participating nuclei in the reaction o)m'^ and superscripts indicate the specific 
excited states. Experiments measure <^i^ijj o; Eij), summed over all excited states of 
the final nucleus, with the target in the ground state. Target states /i in an astrophysical 
plasma are thermally populated and the astrophysical cross section cr,*(j, o) is given by 

E^(2Jf + l)exp(-£;fAT) ■ 
The summation over u replaces T^{E^ J, tt) in Eq. |l| by the total transmission coefficient 
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Here Sm,o is the channel separation energy, and the summation over excited states above 
the highest experimentally known state is changed to an integration over the level 
density p. The summation over target states /i in Eq. has to be generalized accordingly. 

In addition to the ingredients required for Eq. W like the transmission coefficients 
for particles and photons, the width fluctuation corrections W{j, o, J, tt) have to be 
employed. The important ingredients of statistical model calculations as indicated in 
Eqs. through ^ are the particle and 7-transmission coefficients T and the level density 
of excited states p. Therefore, the reliability of such calculations is determined by the 
accuracy with which these components can be evaluated. 



3. Level Densities 

3.1. The Back-Shifted Fermi-Gas Model 

Considerable effort has been put into the improvement of the input for the statistical 
Hauser-Feshbach calculations (e.g. Cowan, Thielemann & Truran 1991). However, the 
nuclear level density has given rise to the largest uncertainties in the description of 
nuclear reactions. For calculating the level densities in that context one does not only 
have to find reliable methods, but also computationally feasible ones. In dealing with 
thousands of nuclei one has to resort to simple models in order to minimize computer 
time. 

Such a simple model is the non-interacting Fermi-gas model introduced by Bethe 
(1936). Mostly, the back-shifted Fermi-gas description, assuming an even distribution of 
odd and even parities, is used (Gilbert & Cameron 1965): 

p{U,J,iT) = ^f{U,J)p{U) , (4) 



with 



^^^^^71^ 12aV4 ' nU,J) = ^^e.p^ j (5) 

The spin dependence is determined by the spin cut-off parameter a. Thus, the level 
density is dependent on only two parameters: the level density parameter a and the 
backshift S, which determines the energy of the first excited state. 

Within this framework, the quality of level density predictions depends on the re- 
liability of systematic estimates of a and 5. Gilbert & Cameron (1965) were the first 
to identify an empirical correlation of the level density parameter a with experimental 
shell corrections S{N, Z) and of the backshift 5 with experimental pairing corrections. 

An improved approach has to consider the energy dependence of the shell effects 
which arc known to vanish at high excitation energies (Iljinov ct al. 1992), i.e. the 
thermal damping of shell effects. Although, for astrophysical purposes only energies 
close to the particle separation thresholds have to be considered, an energy dependence 
can lead to a considerable improvement of the global fit. This is especially true for 
strongly bound nuclei close to magic numbers. 

We use the following description (initially proposed by Ignatyuk, Smirenkin & Tishin 
1975; Ignatyuk, Istekov & Smirenkin 1979) for an energy-dependent level density pa- 
rameter a: 

a{U,Z,N) = ~a{A) l + C{Z,N)l^^^ , (6) 

where 

a{A) =aA + 13 A^'^ (7) 



and 

/(;7) = 1 - exp(-7t7) . (8) 

The shape of the function f{U) was found by approximation of numerical microscopic 
calculations based on the shell model. Thus, wc arc left with three open parameters, 
namely a, /3, and 7. The values of these parameters can be determined by fitting exper- 
imental data. 

Important input for these fits arc the so-called shell correction C{Z,N) and the 
pairing gap A. These have to be taken from microscopic or macroscopic-microscopic 
nuclear mass formulae. Instead of assuming constant pairing or a fixed dependence on 
the mass number A. wc directly determine the pairing gap A from mass differences of 
neighboring nuclei. Thus, for the neutron pairing gap An one obtains (Wang et al. 1992) 

An{Z, N) = ^ [E^{Z, N-2)- 3E'^{Z, N - 1) + 3E^{Z, N) - E^{Z, N + 1)] , (9) 

where E^{Z,N) is the mass for the nucleus {Z,N). Similarly, the proton pairing gap 
Ap can be calculated. 



Previous attempts to find a global description of level densities used shell corrections 
derived from comparison of liquid-drop masses with experiment (C = S = M^xp — A-Zld) 
or the "empirical" shell corrections C{Z, N) = S{Z, N) given by Gilbert & Cameron 
(1965). A problem connected with the use of liquid-drop masses arises from the fact 
that there are different liquid-drop model parametrizations available in the literature 
which produce quite different values for S. 

We realized that one has to reconsider the meaning of the parameter C in the level 
density formula (Eq. ^ . Actually, it should describe properties of a nucleus differing from 
the spherical macroscopic energy (i.e. mass) , including properties which are vanishing at 
higher excitation energies. Therefore, the parameter C should rather be identified with 
the so-called "microscopic" correction i?niic • The mass of a nucleus with deformation e 
can then be written as 

Mthco(e) = £^mic(e) + £;,„ac (sphere) . (10) 

Alternatively, one can write 

Mthco(e) = i^mac(e) + Es+pie) , (11) 

with i?s+p being the shell-plus-pairing correction. The confusion about the term "mi- 
croscopic correction" , which is sometimes ambiguously used in literature, is also pointed 
out in MoUer et al. (1995). Thus, the above mentioned ambiguity of different liquid-drop 
parametrizations follows from the inclusion of partially deformation dependent effects 
into the macroscopic part of the mass formula. 



3.2. Results 

In our study we utilized the most recent mass formula by MoUer et al. (1995), consis- 
tently also taking the microscopic corrections of that mass formula (calculated in a folded 
Yukawa shell model with Lipkin-Nogami pairing) in order to determine the parameter 
C{Z, N)=E^ic- The backshift S was calculated by setting <5(Z, Af)=A„(Z, iV)-|-Ap(Z, N) 
and using Eq. ^. In order to obtain the parameters a, /3, and 7, we performed a fit to 
experimental data on s-wave neutron resonance spacing of 278 nuclei at the neutron 
separation energy. The data were taken from the compilation by Iljinov et al. (1992). 

As a quantitative overall estimate of the agreement between calculation and experi- 
ment, one usually quotes the averaged ratio 
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with n being the number of nuclei for which level densities p are experimentally known. 

As a best fit (when including known experimental masses where known), we obtain 
an averaged ratio g = 1.5 with the parameter values a = 0.1336, f3 — —0.06712, 7 = 
0.04862. The ratios of experimental to predicted level densities for the nuclei considered 
is shown in Fig. |l|. As can be seen, for the majority of nuclei the absolute deviation is 
less than a factor of 2. This is a quite satisfactory improvement over the theoretical level 



densities used previously. Such a direct comparison was rarely shown in earlier work, 
mostly the level density parameter a, entering exponentially into the level density, was 
displayed. Closely examining the nuclei with the largest deviations in our fit, we were 
not able to find any remaining correlation of the deviation with separation energy (i.e. 
excitation energy) or spin. Different combinations of masses and microscopic corrections 
from other models (droplet model by Myers and Swiatecki, Cameron-Elkin mass formula 
and shell corrections, experimental values) were also tried but did not lead to better 
results. 

As the parameters were fit to experimentally accessible nuclei close to the drip line, 
the question arises how reliable such a level density description can be for astrophysically 
interesting nuclei far from stability. Since the microscopic corrections (and the masses) 
are taken from a nuclear-structure model this leads to the question of the reliability of 
that model far off stability. There is, of course, no final conclusion regarding that point 
yet. However, recent investigations in astrophysics and nuclear physics have shown the 
robustness of the MoUer et al. (1995) approach (Thielemann et al. 1994). Recently 
improved purely microscopic models have exhibited similar behavior towards the drip 
lines (Dobaczewski et al. 1994), but there are no large scale calculations over the whole 
chart of nuclei available yet. Therefore, the model used here is among the most reliable 
ones available at present. 

4. Applicability of the Statistical Model 

With the help of a level density description one is also able to make statements on when 
and where the statistical model approach is valid. Generally speaking, in order to apply 
the model correctly a sufficiently large number of levels is needed in the relevant energy 
range which can act as doorway states to form a compound nucleus. In the following we 
will discuss this with the aid of the above level density approach for neutron-, proton- 
and a-induced reactions. This section is intended to be a guide to a meaningful and 
correct application of the statistical model. 

4.1. Neutron-Induced Reactions 

In the case of neutron-induced reactions a criterion for the applicability can directly 
be derived from the level density. For astrophysical purposes the projectile energies are 
quite low by nuclear physics standards. Therefore, the relevant energies will lie very close 
to the neutron separation energy. Thus, one just has to consider the level density at this 
energy. As a rule-of-thumb it is usually said that there should at least be 10 levels per 
MeV for reliable statistical model calculations. The level densities at the appropriate 
neutron separation energies are shown in Fig. ^ (note that therefore the level density is 
plotted at a different energy for each nucleus). One can easily identify the magic neutron 
numbers by the drop in level density, as well as odd-even staggering effects. A general 
sharp drop is also found for nuclei close to the neutron drip line. For nuclei with such 
low level densities the statistical model cross sections will become very small and other 
processes might become important, such as direct reactions. 
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Fig. 1. Ratio of predicted to experimental (Iljinov et al. 1992) level densities at the neutron 
separation energy. 
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4.2. Charged-Particle Induced Reactions 



Contrary to neutron-induced reactions, the knowledge of the level density at the particle 
separation energy is not sufficient to say something about the reliability of the statistical 
model because in this case also the Coulomb barrier has to be taken into account. 
By folding the Maxwell-Boltzmann velocity distribution of the projectiles at a given 
temperature with the penetrability through the Coulomb barrier, the so-called Gamow 
peak can be derived, in which most of the reactions will take place. Location and width of 
the Gamow peak depend on the charges of projectile and target, and on the temperature. 

Similar to the rule-of-thumb given above for neutron-induced reactions, an estimate 
can be made of how many resonances have to be inside the Gamow peak in order to still 
be able to numerically solve the integral when calculating the reaction rate. Comparison 
with experimental cross sections gives the result that one can assume a value of about 3 
resonances (levels). However, there are comparisons to experimental data (Van Wormer 
et al. 1994) which have shown that in certain cases the experimental values could be 
reproduced with even fewer levels. Nevertheless, to be on the safe side we assumed 3.5 
levels within the width (FWHM) of the Gamow peak as the necessary minimum in the 
plots shown in Figs. || and ^. 

For proton-induced reactions Fig. || should be used, whereas Fig. |^ is for a-induced 
reactions. These plots do not directly display level densities but rather show (again in 
a logarithmic scale) the temperatures at which the statistical model can be used. That 
means that for each nucleus (it is always the compound nucleus which is plotted for the 
reaction) the temperature is given beyond which the above described condition of 3.5 
levels is fulfilled. This enables one to directly read from the plot whether the statistical 
model cross section can be "trusted" or whether single resonances or other processes 
(e.g. direct reactions) have also to be considered. (However, this does not necessarily 
mean that the statistical cross section is always negligible in the latter cases, since the 
assumed condition is quite conservative). 

5. Summary 

We were able to improve considerably the prediction of nuclear level densities by employ- 
ing an energy-dependent description for the level density parameter a and by correctly 
including microscopic corrections. All nuclei can now be described with a single parame- 
ter set consisting of just three parameters. The globally averaged deviation of prediction 
from experiment of about 1.5 translates into somewhat smaller deviations of the final 
cross sections due to the dominance of transitions to states with low excitation energies. 
This will also make it worthwile to recalculate the cross sections and thermonuclear 
rates for many astrophysically important reactions in the intermediate and heavy mass 
region. 

We also presented a "map" as a guide for the application of the statistical model for 
neutron-, proton- and a-induced reactions. The above plots can give hints on when it 
is safe to use the statistical model approach and which nuclei have to be treated with 
special attention at a given temperature. Thus, information on which nuclei might be 
of special interest for an experimental investigation may also be extracted. It should be 
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noted that wc used very conservative assumptions in deriving the above criteria for the 
apphcabihty of the statistical modeL 
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